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Abstract 

Within the framework of two-sublattice Mitsui model with taking into account the shear strain 
£4 and the diagonal strains S2 and £3, a dynamic dielectric response of Rochelle salt X-cuts is 
considered. Experimentally observed phenomena of crystal clamping by high frequency electric 
field, piezoelectric resonance and microwave dispersion are described. It is shown that the lowest 
resonant frequency is always associated with the £4 shear mode 
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1. Introduction 

Crystals of Rochelle salt have been attracting an interest of physicists due to their practical 
applications in past, and now from the fundamental point of view mostly. In contrast to most of 
the known ferroelectrics, in Rochelle salt the ferroelectric phase exists only in a temperature in- 
terval between two second order phase transitions at 255 and 297 K. Spontaneous polarization Pi 
is accompanied by shear strain s\. The ferroelectric phase is monoclinic (P2i 11); both paraelec- 
tric phases are orthorhombic (P2i2i2i); all phases are non-centrosymmetric and piezoelectric. 

Dynamical dielectric response of Rochelle salt exhibits several dispersion regions. Those 
are related to the domain wall motion JTJ, [3, S] (below 1 kHz), piezoelectric resonance 0, Si 
(between 10 kHz and 10 MHz), microwave relaxation |7[], and the submillimeter (100-700 GHz) 
resonances JH]. While the very low and very high frequency dispersions are relatively well 
studied, the presence of the piezoelectric resonance dispersion in Rochelle salt is acknowledged 
at best. Quantitative data is available mostly for the temperature dependence of the first resonance 
frequency whereas very little information has been obtained |4] about details of the 

temperature or frequency dependence of the dielectric permittivity, mode assignment for different 
crystal cuts, etc. 

Behavior of Rochelle salt is usually described within a two-sublattice Ising model with an 
asymmetric double-well potential (Mitsui model 1 110 . The pseudospin dynamics is con- 



sidered within the Bloch equations or Glauber approach methods. 

Rochelle salt is a perfect example of a system, whose dynamic dielectric response cannot 
be correctly described without taking into account the deformational effects. Their influence is 
revealed, in particular, in the phenomena of piezoelectric resonance and crystal clamping by a 
high-frequency measuring field, none of which can be obtained within theories based on under- 
formable versions of the Mitsui model iTToL 121. Such theories yield a diverging relaxation time 
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at the Curie point and, as a result, incorrect temperature behavior of the microwave permittivity 
near the phase transitions. 

In II 1311 the dynamic dielectric permittivity of Rochelle salt has been calculated, using the 
model with piezoelectric coupling 11411 with the shear strain £4, for the entire frequency range 
from the static limit (in the ferroelectric phase from about 1 kHz) to 10 11 Hz, including the 
piezoelectric resonance region. For a coupled dynamics of the shear strain £4 - pseudospin 

I — L 

system, the standard methods of description of the lattice strain dynamics [5] based on Newtonian 
equations of motion has been combined with the Glauber approach to pseudospin dynamics. 
Evolution of the dielectric permittivity from the static free crystal value via the piezoelectric 
resonances to the clamped crystal value and to the microwave relaxation has been described. 
Recently, in II 1 511 it has been pointed out that boundary conditions in 111 311 were not set correctly, 
which resulted in the underestimated values of the resonant frequencies; a correct equation for 
the resonant frequencies related to the shear £4 mode has been obtained II 1 511 . 

In the paraelectric phases in Rochelle salt the longitudinal field E\ excites only the shear 
mode £4, as du is there the only non-zero piezoelectric coefficient associated with E\. However, 
in the ferroelectric phase it can also excite the extensional modes associated with the diagonal 
strains s\, £2, and £3 via the non-zero coefficients d\\, dn, dy$. The contributions of the ex- 
tensional modes to the dielectric permittivity of Rochelle salt at frequencies far from the piezo- 
electric resonance region are not expected to be crucial, due to smallness of du (i = 1,2,3) in 
comparison with du. The presence of these modes, however, changes the permittivity in the 
resonance region, at least giving rise to additional resonance peaks, and this should be explored. 

In this paper we follow the same approach that has been used previously 01311 for the case 
with only one shear mode. The modification of the Mitsui model of 01611 that takes into account 
both the shear strain £4 and the diagonal strains is exploited. The drawbacks of the previous 
calculations related to setting the boundary conditions are removed. The expression for the 
dynamic dielectric permittivity and equations for the resonant frequencies of Rochelle salt X- 
cuts is obtained. 



2. System thermodynamics in presence of diagonal strains 

For the sake of the reader's convenience, we shall present here the expressions for the related 
to the shear strain £4 thermodynamic and physical characteristics of Rochelle salt obtained within 
the modified two-sublattice Mitsui model with the shear strain £4 and with the diagonal strains 

M. 

The system behavior is described in terms of the following linear combinations of the mean 
values of the pseudospins belonging to different sublattices 

£ = ^«0-?l> + <0>y2», CT = i«cr 9l ) - (o- 9 2»; 

£ is the parameter of ferroelectric ordering in the system. 

The thermodynamic potential of the model 01611 within the mean field approximation reads 

, J + K 2 J-K 2 21n2 1, , y + S , y-6 
giE{T) = U seed + + — ^— o~ -p p In cosh — ^— cosh — ^— , (1) 



where U see d is the phenomenological part of the thermodynamical potential, representing the 
energy of the host lattice of heavy ions which forms the asymmetric double-well potentials for 
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the pseudospins (see /3 = 1/ksT, £g is the Boltzmann constant, and 

(J + K \ IJ-K \ 

y=/3\——t-2it, i s 4 +fi l E 1 ), (5 = — — cr + Aj. 

Here n\ is the effective dipole moment. The model parameter 1//4 describes the internal field 
created by the piezoelectric coupling with £4. It is assumed that a longitudinal electric field E\ is 
applied. 

The parameters J, K are the Fourier-transforms (at k = 0) of the constants of interaction 
between pseudospins belonging to the same and to different sublattices, respectively. They, along 
with the double well potential asymmetry parameter A, are assumed 10] to be linear functions 
of the diagonal strains 

3 3 
J ± K = J ± K Q + 2 J] ifisi, A = A + J] <A3,£/- (2) 

i=l i=l 

The stress-strain relations and polarization have been obtained in the following form 

7=1 

cr 4 = c^e 4 - e\ A E x + 2^, (3) 
v 

P 1 = e° u s 4 +xf 1 E 1 +^. (4) 

cr, are the components of the stress tensor in Voigt notations; cf? are the "seed" elastic constants; 
a 1 ] are the "seed" thermal expansion coefficients. 



3. Vibrations of X-cuts of Rochelle salt 

We consider vibrations of a thin rectangular L y x L z plate of a Rochelle salt crystal cut in the 
(100) plane (X-cut) induced by time-dependent electric field E\ = E\ t exp(iajt). This field gives 
rise to the shear strain £4 at all temperatures, as well as to the diagonal strains s\, £2, £3 in the 
ferroelectric phase. We take into account the in-plane vibrational modes, allowed by the system 
symmetry, and neglect the out-of-plane mode associated with e\ . 

Dynamics of the pseudospin subsystem will be described within the Glauber approach [ 17], 
where the kinetic equations for the time-dependent averages £ and cr have the form 11411 

- a— J = t- i[tanh -(y + 6) + tanh -{y - 6)}, 
at 1 1 1 

d 1 1 1 

- a ^i°~ = o [ tanh o(r + - tanh ~(y-o)l (5) 

at 2 2 2 

Here a is the parameter setting the scale of the dynamic processes in the pseudospin subsystem. 
The form of (0 is not affected by inclusion of the diagonal strains into consideration. 

Dynamics of the strains will be described by the standard method, using classical (Newto- 
nian) equations of motion 1 18] of an elementary volume 



where p — 1 .767 g/cm 3 is the crystal density, 77, are the displacements of an elementary volume 
along the axis x,, a,u are components of the the mechanical stress tensor. Relevant to our case 
are the displacements 772 and 773, giving the strains 

_ dijz _ drfr _ dr]2 + d/73 

cty ' " 3z ' dz 

At small deviations from the equilibrium the dynamic variables cr, and sj (or 77,) can 
be presented as sums of the equilibrium values and of the fiuctuational deviations, while the 
deviations are taken to be in the form of harmonic waves 

f = I + 6(y, z) exp(;wf), cr = a + a,(y, z) exp(zW), 77,- = fj t + rj it (y, z) exp(iwf)- 

Equations ©, and © can be expanded in terms of these deviations up to the linear terms. 

For £ and a we obtain the same equations that follow from the condition of the thermodynamic 
potential (Q} extremum lfl6ll . 

From (O-© we get the following constitutive equations 

3 1-1 1 

o-i(y,z) = V cf°e jt (y,z) - -tfrf&fy, z) - -ifrjaa,(y,z) - -i/f3iO~t(y,z), 

i—i ' V V V 

7=1 

(T 4 (y, z) = c™e 4t (y, z) - e° u E lt + 2^ft(y, z), (7) 
J>i t (y,z) = e° u e 4t (y,z) +XuE lt + ^? t (y,z). (8) 

(for i = 1 — 3). Taking these into into account, from (0, and © we get the following system of 
equations 

d 

-a—l; t + d\%, + a 2 cr t + a n E2t + fl33£3r + 034^4; = OwMlt, 

at 

d 

-a—a, + big, + b 2 a, + b^s lt + b^si, + b 34 e 4 , = b 02 E u , 
dt 

d 2 T] 2t _ EO&mt , js&mt EO,5 2 772, #W 

P dfi - 22 df 23 8y8z 44 ( dz 2 dydz } 

2i/> 4 <9£ l r , + £ <9£ da, 

v oz v ay ay 

J 2 m, =4 o^L + + ,0(^73, + fl^ + (9) 



df 23 3y3z 33 ^2 -M4V 5y2 - dydz ' 

+ T [^T" + (<A^CT + *A33)^— 

v cry v - dy oz 



with 



J + K K-J 

a\ = -1 +/3 /ti, a 2 = P 4 /t 2 , 

032,3 = ^[Pl ( / , 2,3 _/i 2(i/'2,3 ,5 " + , /'32,3)], ^34 = -ffyiM, 0®i = --^A\ 

J + K J-K 
b!=-fi—j-A 2 , b 2 = -l+/3——A 2 , 

^32,3 = -2^2^2,3 ~ ^l^^ 5 " + ^32,3)L b 34 = /3(/> 4 ^2, ^02 = ^y"^2- 
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(Hereafter it is implied that the deviations are functions of y and z.) We introduced the following 
notations 

A\ — 1 - £ 2 - cr 2 , A 2 — 2^cr. 

Solving the two first equations of (0 with respect to g, at 772^ = i]3t = (regime of a mechani- 
cally clamped crystal), substituting the result into ([SJ, and differentiating that with respect to the 
field, we find the dynamic dielectric permittivity of a clamped crystal 

Ru 2 

e E u (cj) = e$ + An^F^aai), sf, = 1 + 4^?. (10) 



where 



^ 1W = — — -j— — - — - — , <p 1 =2-—A u 
(iaa>Y + {iau>)<pi + <p 2 2 

<P2 = l- —M -p — — — (Aj - A 2 ), if 3 = Ai + (A l - A 2 ) 

1 to 4 



This is the same expression that has been obtained previously 11311 for the Mitsui model with the 
shear strain £4, but without the diagonal strains. 

In the regime of a mechanically free crystal, the two first equations of (01 give 

&(y.z) = -^-Fi(aa))Eu - P^AF\{auj)E A ,{y,z) - 



-F 42 (aaj)e 2t (y,z) - -F^ 3 (aw)e 3 ,(y,z), 
2 



o~t(y,z) = —F%(aaj)E lt -pip 4 F'[(au))s4t(y,z) 



~ ^F% 2 (au)s2,(y, z) - -F% 3 (au)s 3l (y, z), (11) 



where 



FfXaoj) 



tf 4j + (zW)[f/li^J - (^<x + i// 3 j)A 2 ] 



4 J y D{cS) 



D(co) 

A 2 (l + iato) 



D(a>) = (iaco) 2 + {iau)tp\ + ip 2 , 

<P4i = (A, + ^3 - (if/Jv + ^3i) h., <P5 = M - ^—j—(A\ - A\). 



Hence, the field-dependent part of polarization is 



Pu(y,z) = 

where 



xfi + ^F^ato) 



Ei t + 2^ e ii(^)si,(y,z), (12) 



i=2,3,4 



ei2,3(w) = — Fl 23 (aco), eu(u>) = e u — F\(au) (13) 



are dynamic piezoelectric coefficients. 

The observable dynamic dielectric permittivity at constant stress xli ( w ) is expressed via the 
derivative from the polarization averaged over the sample volume 

£n(oj)=l + —— dy dzPify,z) = 
if Jo Jo 



L y L z dE x 

dE u 



s s n (o))+4n J] e u (co)—e it (y,z), (14) 

1=2,3,4 ' 



it = -r~r \ d y f dze it (y,z). 

L y L z Jo Jo 



where 

<°, 

The remaining problem is to find Si t (y, z). 

Taking into account ( fTTT ). from the two last equation of (O it follows that 

-pa> m< = c 22(^)-^r + 2c 2^Qyfc + c 44( w )-^r + 

2 ~E , ^d 2j l2t r~E , -> ~E , m ^Tfc ~E , ^^ Z T] 2 , 

-pOJ T]3, = C 24 (w)-^2- + [C 23 (W) + C 44 (W)]^- + C 34 (w)-^- 2 - + 

~e, ,d 2 m, n .E, ,d 2 mr _ E( a 2 m , 
+ <%M-£r + 2c 3 4 M^ + 4(^)^r- (15) 

which is nothing but the Christoffel equations, with the frequency-dependent elastic constants 
cfj(a>) given by the model expressions 

-E, , ^4 rf , . 

c, 4 (w) = —F^iaco), 

~E,. ,\ £0 P 



c* 4 (cj) = eg - ^Fi(aw). (16) 

Their frequency variation is perceptible only in the region of the microwave dispersion of the 
dielectric susceptibility. However, in the piezoelectric resonance region, which is expected to be 
in the 10 4 - 10 7 Hz range, depending on temperature and sample dimensions, well separated from 
the region of the microwave relaxation, cfj{a>) and e,y(w) < TT~3T > are very close to the corresponding 
static quantities (see lfl6l[l9ll ) and coincide with them at co — > 0. 

We can rewrite $15[ in a different form, where it will be more convenient to set the boundary 
conditions. We differentiate (fT~5b with respect to y and z, transforming it into three equations for 
the strains si t , S3,, S4, instead of the displacements j] 2t and 773, 



2 -£ , . d 2 s 2 , , E d 2 s 2 , _ E . . d 2 e 2t 
j e 2 , = c 22 ( W )— + c 2A (oS)— + c 44 (co) — 



, E .d 2 S 4t r _£. . „ E , .,d 2 S3, _ E .d 2 E3 t 



2 ~E , ,d E2t r-E / \ ~E / M 5 e 2f .<3 2 £ 4 / 

-pw-e 3 / = cf 4 (w)^-^ + [cf 3 (w) + c^a/)]-^- + 4 4 (to)-^- + 

-e , s <3 2e 3f _g 3 2 e 3f _ £ , . d 2 s it 
+ c AA {u» w+ c M (u»— + c^u)—, 

-pa> z s 4 , = c^iu)-^- + [c£ 2 (w) + c m( w )]^; + [2cf 4 (w) + ^(w)]-^- + 

+ [<&(£<>) + 2c f 4 ( w )]^l- + t c 23( w ) + C 33( w )]t3T + c 34( w ) 



The boundary conditions for e„ follow from the assumption that the crystal is simply sup- 
ported, that is, it is traction free at its edges (at y = 0, y — L y , z = 0, z = L z , to be denoted as 

0"ils = <r 2 \x = o- 3 | s = <r 4 k: = 0. (18) 



In our previous consideration 1 13] this condition was fulfilled at the corners of the crystal plate 
only, not at all its edges. Eventually that led to an incorrect expression for the resonance frequen- 
cies. 

Substituting ( fl~8b into the constitutive relations (Q) and using ( TTTb . we obtain the boundary 
conditions for the strains in the following form 

Enk = £;o = du(a))Ey, (19) 

where 

4 

d u (a>) = ^ 4^ e ^' ( 20 > 
/=i 

eii(a>) is given by ( fT3l , and s?(o>) are the elements of a matrix inverse to cfj(u>) determined in 
Let us consider first the case of the paraelectric phases (at £ = 0). Then 

620 = £30 =0, (21) 

as expected from the symmetry considerations. As we shall show later, from this boundary 
condition it follows that 

sit(y,z) = S3t(y,z) = 0, 
at all (y, z). Then the system dTTb reduces to a single equation for S4 t (y, z) 

2 ~e , M-d 2 ^/ d 2 s At 
-poj s 4 , = c 44 (a))[-^- + (22) 



Its solution is 



S 4t (y, Z) = £40 + (23) 
16 to 2 . n(2k+l)y . n(2l + V)z 



sin sin • 



k ^ Q (2k + 1)(2/ + 1 )n 2 (a) W)2 _ U 2 L y L z 

7 



with oJ: ? given by 



p 



(2k + l) 2 (2/+1) 2 



3 



L 2 



(24) 



Since in the paraelectric phases 



ei4(a>) 

£40 = _ F . . E\ t = du(ci))Eu, 
C44M 



we have that 
where 



£4 f = R4{cS)d\<\((jS)E\ t , 
64 

7? 4 (w) = 1 + > : 7— 

4V ; Zj (2Jk + 1) 2 (2/ + 1) 2 7T 



64 w 2 



^ ) (2fe + l) 2 (2/ + l)V (w (4) ) 2_ w 2- 

Thus the dielectric permittivity reads 

£n(w) = £n(w) + 47r/?4(w)ei4(<y)ofi4(w). 



(25) 



(26) 



Let us analyze the above results. In the static limit (d — > 0, R^{oS) — > 1) from (1261 1 we obtain 
the static permittivity of a free crystal (see B19I1 ): in the high frequency limit (2jfc/=o64/[(2A: + 
1) 2 (2Z + l) 2 7r 4 ] = 1, and Ra{(l>) — * 0) we get a dynamic permittivity (TTOb of a mechanically 
clamped crystal, exhibiting relaxational dispersion in the microwave region. Thus, eq. (|26*T > 
explicitly describes the effect of crystal clamping by high-frequency electric field. 

In the intermediate frequency region, it has a resonance dispersion with numerous peaks 
ar frequencies where Re[R4(a>)] — > 00. In this frequency range we can neglect the frequency 
dependence of c^ 4 (u>) and reduce the equation for the resonance frequencies (l24T i to an explicit 
expression by putting in it c£(a>) — > c; 



u 44- 



Comparing d24l > to the expression obtained previously IU3I1 for a square X-cut 



a> k 



K(2k+l) 



we can see that the incorrectly set boundary conditions [13] led to the V2 times smaller lowest 
resonance frequency than the correct one. However, the low and high frequency limits of the 
permittivity [13] (the static value and the clamped values with the relaxational dispersion in the 
microwave region) were correct. 

Now we shall proceed to the case of the ferroelectric phase. The system of second-order 
partial differential equations ([T7l > will be solved numerically using the finite element method. 
However, its main features, such as most of the resonant frequencies, including the lowest one, 



can be obtained analytically, if we take into account the fact that C22, C33, c^, C23 » cf 4 , cf 4 . 
Neglecting in dT7T > the terms proportional to cf 4 , c^ 4 (in the paraelectric phases Cj 4 = cf 4 = 



exactly) we get 



-poj e 2t = c 22 (uj)— -r- + c A4 (u>)— y- + [c 44 («) + c 23 (w)] 



dy 2 



dz 2 



d 2 s 2 , 



-pa> e 3 , = [c 23 (u>) + c^Cai)]-^- + 



c 44 (w) 



d 2 S3, 

dy 2 



dy 2 ' 
d 2 s 3l 
dz 2 ' 



9 F F ^^^2/ F F 3^ 

-pco l s 4 , = [cgj(G>) + c^(w)]— + [cf 3 (w) + cf 3 (w)]r^- + 



dydz 
]• 



dydz 



(27) 



The system is partially split, with the two first equations depending on and £3, only. We look 
for the solutions in the form of series 



Sit = BiO 



+ > Z). sin sin , i = 2,3, 



4. 



(28) 



It is easy to verify that the boundary conditions (fT9l l are satisfied. 

Substituting ( l28T l into the two first equations of ( l27l i and noting that we can write that 



2 7r 2 (2/fc + 



16 



(2Jt + 1)(2/ + 1) 



. n(2k+ l)y . tt(2/ + l)z 
■ sin sin , 



we obtain 



with 



m = 



16 



7r 2 (2fc+l)(2/+l) A (f ) ( w ) 



-di 3 (w)(cf 3 (w) + c^(w)) 



7,E 



n 2 {2k + l) 2 



16 



, ^ 2 (2fc+l) 2 , g ^ 2 (2/+l) 2 2 

C 44 (w) ^5 + c 33( w ) P w 



L 2 



7r 2 (2fe+l)(2Z+l) A (f ) ( w ) 



-c/i 2 (<y)(cf 3 (w) + c^(w)) 

T 2 



n 2 (2l+\) 2 



'_ E n 2 (2k+l) 2 _ E 7i 2 {2l+\) 2 



c 22 (w)- 



1 



+ C44M- 



L 2 



L 2 

2 \ 
— pco 



cm f n 2 {2k+\f F n 2 (2l+\) 2 , 

= [£f 2 M L2 + L2 "PW 2 ] X 



x[c^ 4 (w) 



7r 2 (2fc + l) 2 



+ cf 3 (w) 



7T 2 (2/+ l) 2 



z 2 ' • !! "" /,: 



-pa>]- 



,~e, s^~e, „ 2 n 2 (2k+l) 2 n 2 (2l+iy 
-(c 23 (w) + c 44 (w)) 



3 



(29) 



2 



(30) 



0. 



Note that owing to the boundary conditions (t2Tb . in the paraelectric phases D 
It means that the extensional modes associated with the strains £2 and £3 are not excited by the 
longitudinal field E\, which is expected from the symmetry considerations for the orthorhombic 
point group of Rochelle salt. The same results are obtained by the numerical finite element 
method calculations for the complete system (fT7l ). 

Substituting the found S2t(y,Z) and £3,(y, z) from ( |28| ) with ( 1291 into the third equation of 
(TZTT i and then reexpanding two first terms in its right-hand-side in series over sin —j^- sin 2^2, we 
obtain that 



Df = 



l6du(tL>) 



n 2 (2k+ \)(2l+ 1) 



TT) + 2^ e » 



ipkmpln 



OJ 



(4?) 2 -w 2 ' 



(31) 



with 



Q mn = [cUu>) + c% 3 (")]D™ + [g* (w) + («j)]D™, (32) 

and 



1 

Pkm - 



m + k+\ 



if k and m are of the same parity, and 



Pkn 



otherwise. 

The averaged over the sample volume strains occurring in the expression for the permittivity 
are then equal 

s~t = Ri{o))d\i{cS)Ei, (33) 

with 

Ri(co) = 1 + > -r '■ (34) 

n 2 {2k+ 1)(2/+ l)d u (to) 

and Df given by d29| i and (fJTJi. Finally, the permittivity is 

s^ioS) = efj(w) +4n eiKw)rfi;(w)/?;(a/). (35) 

'=2,3,4 

We are in a position now to determine the resonant frequencies of the permittivity in the 
ferroelectric phase. Those occur at Z?,(o>) — > °o and are of two types. The first type of resonances 
is given by /?2,3(w) — > oo, that is by equation 

Ag 3) (w) = 0, (36) 

and it is associated with the extensional modes of £2 and £3. Solutions of Q61> . to be denoted as 

(23) 

co), , exist at all temperatures, but the corresponding modes are not excited in the paraelectric 
phases. Therefore, these resonances are present in the ferroelectric phase only. On the other 
hand, R^u) — > 00 both at <±> — w^ 3) and at a> — or,) given by (1241 1. The resonances given by 
(l24l originate from the shear vibrational mode associated with £4 and persist in the paraelec- 
tric phases. Such a division is, however, artificial, as the modes are coupled. All three strains 
calculated numerically from the complete system (fTTT i have resonances at the same frequencies. 



4. Numerical analysis 

The set of the model parameters, providing a fair description of dielectric, piezoelectric, and 
elastic characteristics of Rochelle salt, its microwave permittivity, as well as thermal expansion 
of the crystal and the effects of hydrostatic and uniaxial pressures has been obtained in 1 16, 1^1 . 
No additional theory parameters need to be determined apart from those. However, the sample 
dimensions should be specified. In the present paper we shall use L y = 1 .60 cm, L z = 2.45 cm of 
the Rochelle salt X-cut sample, for which experimental data on the resonant frequencies and the 
dielectric permittivity in the resonance region are available [4]. The static (equilibrium) values 

10 



of the dynamic variables f , cr, e, are calculated by minimization of the thermodynamic potential 
(Q} with respect to maximization with respect to cr, and from equations (f3]). Eqs. ( TTTb for the 
strains were solved numerically with the finite element method package FreeFem++ 12011 . The 
solutions were used to find and, hence, the permittivity. 

In figures [1] and [2] we show the spatial distribution of the fluctuational parts of the strains 
Sit(y,z)/Ei, at different frequencies, obtained by solving the full system ( TTTb with the boundary 
conditions ( fT9b . Note that the gray scales are different for each graph. 

The distributions have a single extremum at the sample center at low frequencies (up to the 
frequency of the first resonance); then the extrema multiply. Above the resonances, e,- r (y, z) are 
zeros in most of the plate, only going to the boundary values given by ( fT9l ) within very narrow 
strips near the sample edges. It illustrates the effect of crystal clamping by a high-frequency 
electric field. 

Figure[3]shows the frequency dependence of dynamic permittivity of the Rochelle salt X-cut 
(with L v = 1.60 cm, L z = 2.45 cm) within the entire frequency range of the current model ap- 
plicability. That range does not include the region of the domain-related dispersion below 1 kHz 
Ol, IjXI or the submillimeter (100-700 GHz) region of resonant dispersion [8[]. The experimental 
data shown by open symbols are for the frequencies outside the piezoelectric resonance regions 
of the samples used in the measurements, so the dimensions of those samples are irrelevant. The 
obtained evolution of the permittivity is analogous to the experimental [3] and to the previously 
obtained theoretical lfl3ill5Tl ones: from the static permittivity of a free crystal at low frequencies, 
via the piezoelectric resonance region (10 4 -r 10 7 Hz for this sample dimensions), to the clamped 
crystal value, and, eventually, to a relaxational dispersion in the microwave region. A fairly good 
agreement with experiment is obtained. 

Now we take a closer look at the resonance region. Ability of the simplified system ( 1271 ) to 
describe the resonant behavior of lattice strains and, henceforth, the dynamic dielectric permit- 
tivity of Rochelle salt X-cuts is demonstrated in fig. |H showing the frequency dependence of 
the permittivity in the lower part of the resonant region. One can see that most of the resonant 
frequencies of the permittivity, calculated numerically, using the complete set of equations ( fTTI ). 
are very well reproduced by the resonant frequencies of the simplified system ( 1271 ). 

The temperature dependence of a few lowest resonance frequencies of the simplified system 
(l27l i for this particular X-cut is shown in figure|5] It is seen that at all temperatures the lowest is 
the resonance frequency vf? = oJy* /2n associated with the shear mode at k — I — 0. The lowest 

(23) (23) 

frequency v), — or,, /2tt of the extensional modes (at k = I = 0), active in the ferroelectric 
phase only, is higher at any temperature and at any sample dimensions and always remains finite. 
The shear mode frequencies vi, , on the other hand, go to zero at the Curie temperatures at all k 
and /. 

The agreement with experiment for the lowest resonant frequency is quite good in the fer- 
roelectric phase and gets worse in the upper paraelectric phase. It is, apparently, caused by a 
similar misfit for the elastic constant cf 4 (see Eq. d24"l) and B19I1 ). 

This disagreement is well revealed in the temperature variation of the dynamic dielectric 
permittivity of the Rochelle salt X-cut at 60 kHz shown in fig. [6] (left). There are two distinct and 
well resolved resonance peaks below and above the Curie point. With increasing frequency, these 
peaks move away from the Curie point. They are associated with the £4 shear mode and given 
by vS?. In the upper paraelectric phase, the disagreement between theory and experiment for v res 
seen in fig. [5] resulted in a more than 2 K difference between the theoretical and experimental 
temperatures of the resonant peak. 
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What has not been observed experimentally J3l is that apart from the two discussed peaks, 
in the close vicinity of the transition temperature, the permittivity also has a multitude of other 
resonances, which are higher order resonances of the £4 mode (vf) at k + I > 0). 

Actually, it follows from fig. [5]that at frequencies below the lowest resonant frequency of the 

(23) 

extensional mode (7c) at the Curie points (the threshold frequency which can be found from 
(l36l l at c^ 4 — > 0; it equals 80 kHz for this particular X-cut), all resonances are associated with the 

(23) 

shear £4 mode. Above vL (7c), the resonances associated with the extensional modes appear in 
the ferroelectric phase. 

The temperature curve of the permittivity in a wider temperature range shows (fig. [6] right) 
that the resonances form two packs around each Curie temperature, with the density of reso- 
nances increasing at approaching the Curie points. It can be shown that the pack widths increase 
with increasing frequency, which is caused by the shape of the v?) AT) curves (see fig. O and 

(23) 

by appearance of the v., resonances. Eventually, at some sufficiently high frequency the packs 
overlap. 

5. Concluding remarks 

Vibrations of X-cuts of Rochelle salt crystals and their influence on the dynamic dielectric 
permittivity are analyzed using the modified Mitsui model that takes into account the shear strain 
84 and the diagonal strains e\ , £2, £3 11611 . The system dynamics is described within the frequency 
range, starting from 1 kHz (above the dispersion associated with the domain wall motion) via the 
piezoelectric resonance region and the microwave relaxational dispersion up to about 10 11 Hz. 

Special attention is paid to the piezoelectric resonance region. Explicit expressions for the 
resonant frequencies, associated with the shear mode of £4 and with the extensional in-plane 
modes of £2, £3, of such cuts are derived. They are obtained at neglecting the out-of-plane mode 
associated with s\, as well as the elastic constants c^ 4 and cf 4 . The temperature behavior of the 
resonant frequencies is analyzed; it is shown that the lowest resonance is associated with the 
shear mode at all temperatures. 

The changes in the calculated spatial distributions of the strains with increasing frequency 
visualize the effect of crystal clamping by the high-frequency electric field. Both the shear mode 
and the extensional modes are suppressed. 

It is shown that the resonances associated with the extensional modes appear above a certain 
threshold frequency, and in the ferroelectric phase only, which is consistent with the symme- 
try considerations. In the close vicinities of the transition temperatures, the permittivity has a 
multitude of overlapping peaks, which are higher order resonances of the £4 mode. 
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Figure 3 : Frequency dependence of dielectric permittivity of a X-cut of Rochelle salt at 293 K. Symbols are experimental 
points taken from @l - 01 - V, l2lll - o. The line: a theory. The line and ■ are for Ly = 1.60 cm, L z = 2.45 cm. 
Other symbols correspond to samples of different sizes. 




Figure 4: Frequency dependence of the dynamic dielectric permittivity of a Rochelle salt X-cut at 275 K (left) and 293 
K (right). The solid line is calculated with £/7 found from <17t . A and ■ are the resonant frequencies Vj) and y^ 3 ' of the 
simplified system (27). L v = 1.60 cm, L z = 2.45 cm. 
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Figure 5: Temperature dependence of the lowest resonance frequencies of the Rochelle salt X-cut with Ly = 1.60 cm, 
L- = 2.45 cm. Lines: the theory. Solid lines: v.?; dashed lines: v^ 3 '. Symbols: experimental points of 01. The numbers 
in parentheses are the (k, I) values. 




290 295 300 305 T (K) 24 2 60 280 300 J (K) 



Figure 6: Temperature dependences of dynamic dielectric permittivity of the Rochelle salt X-cut with Z™ = 1.60 cm, 
L z = 2.45 cm at 60 kHz (left) and 800 kHz (right). Solid line: the theory. •: experimental points of [4|]; A are the 
resonant frequencies of the simplified system j27t given by {24). 
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